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Abstract 

Background: Leaf pigment content is an important trait involved in environmental interactions. In order to 
determine its impact on drought tolerance in wheat, we characterized a pale-green durum wheat mutant 
{Tritlcum turgidum L. var. durum) under contrasting water availability conditions. 

Results: The pale-green mutant was investigated by comparing pigment content and gene/protein expression 
profiles to wild-type plants at anthesis. Under well-watered (control) conditions the mutant had lower levels of 
chlorophylls and carotenoids, but higher levels of xanthophyll de-epoxidation compared to wild-type. 
Transcriptomic analysis under control conditions showed that defense genes (encoding e.g. pathogenesis-related 
proteins, peroxidases and chitinases) were upregulated in the mutant, suggesting the presence of mild oxidative 
stress that was compensated without altering the net rate of photosynthesis. Transcriptomic analysis under terminal 
water stress conditions, revealed the modulation of antioxidant enzymes, photosystem components, and enzymes 
representing carbohydrate metabolism and the tricarboxylic acid cycle, indicating that the mutant was exposed to 
greater oxidative stress than the wild-type plants, but had a limited capacity to respond. We also compared the two 
genotypes under irrigated and rain-fed field conditions over three years, finding that the greater oxidative stress 
and corresponding molecular changes in the pale-green mutant were associated to a yield reduction. 

Conclusions: This study provides insight on the effect of pigment content in the molecular response to drought. 
Identified genes differentially expressed under terminal water stress may be valuable for further studies addressing 
drought resistance in wheat. 

Keywords: Chlorophyll content. Oxidative stress. Pale-green mutant, Proteomics, Transcriptomics, Triticum turgidum 
L var. durum, Water stress 



Background 

The leaf chlorophyll content of plants determines their 
capacity to absorb energy from sunlight and is therefore 
considered a measure of photosynthetic potential [1]. 
Plants typically respond to water stress by reducing the 
leaf chlorophyll content, and hence their photosynthetic 
activity [2]. Genotypes that tolerate water stress tend to 
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have more leaf chlorophyll and lower canopy tempera- 
tures when stressed [1]. In contrast, genotypes that are 
susceptible to water stress often have higher canopy 
temperatures [3], and when such plants are well- 
watered, with their stomata fully open, water evapor- 
ation causes leaf cooling [4]. The ability of a genotype 
to remain green under water stress conditions may 
therefore enhance photoassimilation, thus increasing 
productivity and grain filling [5]. In contrast to these 
general statements, a low level of leaf chlorophyll can 
limit the interception of light to avoid oxidative damage 
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caused by excess radiation [6] and could be considered a 
desirable trait under water stress conditions [4]. The 
optimization of photosynthesis in crops under water 
stress therefore requires a productive balance between 
assimilation under favorable conditions and the avoid- 
ance of excess radiation under stress [4]. 

Chlorophylls and carotenoids play an important role 
in photoassimilation and the synthesis of both molecules 
is coupled to chloroplast development [7]. Thylakoid 
organization is coordinated with the expression of nu- 
clear genes, encoding pigment-binding proteins that are 
imported into the chloroplast to assemble the photosyn- 
thetic complexes and the expression of chloroplast 
encoded genes [8,9]. This process is tightly regulated, so 
mutants deficient in chlorophyll content can arise not 
only through the disruption of chlorophyll biosynthesis 
but also via processes related to photosystem assembly 
and light-harvesting activity [10,11]. Such chlorophyll 
mutants have been used to study the function of photo- 
synthetic components as described in rice e.g. [12,13], 
barley e.g. [14-17] and wheat e.g. [18-20]. 

Carotenoids are involved in photoprotection and in 
light harvesting, and different carotenoids play different 
roles in the photosynthetic complexes. For example, only 
P-carotene is present in the reaction centers whereas 
xanthophylls are present in the antennal light harvesting 
complex (LHC) where each protein has a specific carot- 
enoid complement that ensures maximum photosystem 
efficiency [21]. Lutein is involved in light harvesting and 
the quenching of ^Chl states [22] and violaxanthin can 
be replaced by zeaxanthin if the de-epoxidation cycle is 
activated by photosystem saturation following the ab- 
sorption of excess light. Zeaxanthin helps to prevent 
photo-inhibition by dissipating the excess of energy as 
heat, either by direct quenching or by inducing conform- 
ational changes in the LHC proteins [23]. Photo- 
inhibition is an early consequence of excess light and 
other environmental stresses, reducing photosynthetic 
activity and leading to pigment photo-oxidation, the pro- 
duction of reactive oxygen species (ROS) and ultimately 
cell death [24]. Chloroplasts and mitochondria are the 
primary sources of ROS, which damage cells by oxidiz- 
ing nucleic acids, proteins and lipids. However, low con- 
centrations of ROS act as potent signaling molecules 
that coordinate responses to abiotic stress, pests and 
pathogens [25]. The balance between ROS production 
and scavenging by antioxidants can be disrupted by dif- 
ferent forms of stress [26] . 

Wheat is a staple crop that plays a major role in global 
food security. Conventional breeding strategies have been 
used to increase grain yields and stress tolerance, but in 
order to feed a growing population it will be essential to in- 
clude molecular breeding as well as conventional ap- 
proaches to achieve further substantial gains [27,28]. Gene 



expression profQing with microarrays has been used to 
characterize the response of wheat to heat and drought 
stress [29], cold [30], Fusarium infection [31] and phase- 
transition [32], and to investigate the molecular basis of 
grain quality traits [33]. In the context of drought stress, 
Krugman et al. [34] compared resistant and susceptible wOd 
emmer genotypes to investigate the overlapping regulatory 
and signaling processes overrepresented in the drought- 
resistant genotype and to identify candidate genes for 
drought tolerance. Proteomics has also been used to inves- 
tigate the response of wheat to nitrogen [35], salinity stress 
[36,37], low temperatures [38] and drought [39,40]. Poor 
correlation between transcriptomic and proteomic data has 
been reported in some stress experiments although such 
complementary information highlights the multilevel regu- 
lation of responses to complex environments [41,42]. 

Most of the experiments discussed above involved geno- 
types with contrasting phenotypes and different genetic 
backgrounds, an approach that allows many genes and 
stress tolerance strategies to be identified, but makes it diffi- 
cult to characterize specific genotype-phenotype relation- 
ships. This can be achieved in experiments that involve 
genotypes differing at a single locus e.g. [30]. In order to 
control for genetic background effects, we compared a 
pale-green mutant (M) and its wild-type mother variety 
(WT) to determine the impact of leaf pigment content on 
gene and protein expression. We studied the chlorophyll 
and carotenoid content as well as the transcriptomic and 
proteomic profiles of the mutant and wild-type plants 
under well-irrigated (control) conditions and terminal water 
stress. We also tested the performance of the mutant and 
wUd-type plants under irrigated (control) and rain-fed field 
conditions in order to assess the effect of the mutation on 
biomass accumulation and grain yield. 

Results 

Chloroplast ultrastructure 

Wild-type and mutant leaf sections at the tillering stage 
were compared by transmission electron microscopy 
(TEM). As shown in Figure 1, lens-shaped chloroplasts 
were present in the wild-type mesophyll cells whereas 
those in the mutant plants were typically cylindrical. 
The thylakoid grana of the wild-type chloroplasts ap- 
peared to be connected by continuous intergranal lamel- 
lae parallel to the long axis of the chloroplast, whereas 
the mutant grana appeared to be disconnected and dis- 
organized within the chloroplast stroma. Osmiophilic 
lipid globules (plastoglobules) were more abundant in 
the mutant chloroplasts (shown by arrows in Figure IC 
and Figure ID). 

Physiological responses and field measures 

Water stress was applied to plants in the glasshouse by re- 
ducing the soil water content at anthesis. This increased 
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Figure 1 Chloroplast ultrastructure. Transmission electron microscopy (TEM) images of a mesophyll cell from (A) wild-type and (B) pale-green 
mutant leaf (Bar: 5 |jm). (C) Detail of wild-type and (D) mutant chloroplasts, the latter showing a disorganized pattern of unconnected grana and 
several plastoglobules (arrows; Bar: 0.23 nm). 



the negative value of the midday flag leaf water potential 
Q¥w): reduced the leaf conductance to water vapor values 
(gs) and reduced the mean net photosynthetic rate (An) but 
differences between genotypes were not statistically signifi- 
cant (Table 1). In the field, differences in aboveground bio- 
mass were significant only under rain-fed conditions, but 
the grain yield of the mutant plants was reduced by 10.3% 
under irrigated conditions and by 14.3% under rain-fed 
conditions, compared to wild-type plants in the same envir- 
onment (Table 2). 

Pigment content 

The total carotenoid and chlorophyll levels in the flag 
leaves at anthesis under well-watered (control) condi- 
tions were significantly lower in the mutant than the 
wild-type, but differences in antheraxanthin and zeaxan- 
thin levels were not statistically significant (Table 3). 
Under water stress conditions, only the antheraxanthin 
content differed significantly between wild-type and mu- 
tant leaves, whereas the total carotenoid and chlorophyll 
contents were similar in both genotypes. The xantho- 
phyll de-epoxidation cycle appeared to be induced in the 
mutant plants under control conditions. Whereas under 
terminal water stress, wild-type plants increased the de- 
epoxidation rate, mutant plants did not. Therefore, there 



was no statistically significant difference between the ge- 
notypes under water stress conditions (Table 3). 

Transcriptomic profiling 
Microarray quality analysis 

Microarray data were collected using the Affymetrtx 
Wheat Genome Array and normalized using the RMA 
algorithm. The average background was 5.33, well within 
the recommended levels. The percentage of "Present" 
calls ranged from 39.05% to 42.10% among the 61 k 
probe sets. The quality of the biological replicates was 
determined by calculating values among replicates of 
the same sample and the values ranged between 0.988 
and 0.997 with an average of 0.995. Principal component 
analysis highlighted the main sources of total variance. 
The two components explaining 59.3% and 15.8% of the 
variance represented the differences between the treat- 
ments (control vs. water stress) and between the geno- 
types (pale-green mutant vs. wild-type) as shown in 
Figure 2. 

The mutation effect on gene expression was studied 
comparing the expression of the mutant to wild-type for 
each water treatment and indicated along the manu- 
script with the terms induced or repressed genes. In 
addition, the stress effect was evaluated comparing the 



Table 1 Analysis of variance and mean values for physiological traits under glasshouse conditions 



Booting 



Anthesis 



SWC 

(%) 



(Mpa) 



(mol m'^s"^ 



An 

(^mol m"V^) 



SWC 

(%) 



(Mpa) 



An 

(jimol m"^s"^) 



A) ANOVA F-values 

Effect 

Water treatment 
Genotype 

Water treatment x genotype 

B) Mean values 

Water treatment 

Control 

Water stress 

Control 

Wild-type 

Mutant 

Water stress 

Wild-type 

Mutant 



0.36 n.s. 
0.02 n.s. 
0.05 n.s. 



47.1 
48.9 



47.f 



49.0 
48.7 



0.33 n.s. 
0.33 n.s. 
8.22 n.s. 



-0.500 
-0.525 

-0.550 
-0.450 

-0.450 
-0.600 



0.040 n.s. 
0.001 n.s. 
0.100 n.5. 



0.255 
0.255 

0,254 
0,257 

0.259 
0,251 



1.28 n.s. 
1.07 n.s. 
0.97 n.s. 



10.03 
10.49 

10.32 

9.73 

10.50 
10,49 



405 

1.96 n.s. 
0.52 n.5. 



47,2' 
13,5' 

45,4 
48.9 

12,9 
14,0 



218.5 
3.56 n.s. 
0.61 n.s. 



-0.555 " 
-1.205 

-0.580 
-0.530 

-1.250 
-1.160 



90.1 

0.080 n.s. 
0.001 n.s. 



0.368 " 
0.135 

0.365 
0.371 

0.131 
0.139 



23.5 

0.83 n.s. 
2.00 n.s. 



15.94" 

10.85 '° 

14.72 
17.16 

11.11 

1058 



A) F-values of the analysis of variance and, B) means of soil water content (SWC), flag leaf water potential (Tw), stomatal conductance (g^) and net photosynthetic rate (An) for each water treatment and for the wild- 
type and the mutant plants under control and water stress conditions. Water stress started at booting. ***P < 0,001, n,s. P > 0.05. For each variable, water treatment means with different letters are significantly 
different for a Tukey test. 
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Table 2 Yield and biomass under field conditions 



Water regime 


Genotype 


Yield (kg/ha) 


Biomass [g/m^] 


rrigated 


WT 


7913" 


1093" 




Mutant 


7101'' 


1056" 


Rainfed 


m 


5881" 


859" 




Mutant 


5041 


744b 



Average values of grain yield and aboveground biomass were determined at 
anthesis in the wild-type (WT) and the pale-green mutant. Data represent 
means of 3 years of field experiments. Genotype means with different letters 
within a water regime differ at P = 0.05. 



systemic interaction with the environment, abscisic acid 
(ABA) responses, hormonal regulation, stress responses 
and DNA degradation (Table 4), e.g. sequences repre- 
senting RAB (responsive to ABA) proteins, late em- 
bryogenesis abundant protein LEAl and endonucleases 
(Table 5). Probe sets related to cell rescue, defense and 
virulence (functional category 32) were induced or re- 
pressed in the mutant when compared with the wild- 
type. However, none the subcategories was statistically 
significant atP< 0.005. 



gene expression of the water stress to control conditions 
for each genotype and indicated along the manuscript as 
up/downregulated genes. 

Gene expression under control conditions 

Transcriptional profiling at anthesis revealed differences 
in gene expression between the two genotypes under 
control conditions. We identified 23 induced and 15 
repressed probe sets in the mutant (>2-fold difference; 
P < 0.05). The array data for selected probe sets were 
validated by quantitative real-time RT-PCR, showing 
strong agreement between the fold-change values ob- 
served with each method (r = 0.96, P < 0.0001; Additional 
file 1: Figure SI). 

MIPS FunCatDB was used to distribute the genes into 
functional categories, revealing that the induced probe sets 
in the mutant represented chitin catabolism, defense (par- 
ticularly systemic acquired resistance) and cell wall metab- 
olism (Table 4). The most strongly induced probe sets 
represented peroxidase 10 (POXIO), pathogenesis-related 
(PR) proteins including beta-glucanases and chitinases, 
and a TaWRKY45-like transcription factor (Table 5). The 
repressed probe sets represented functions including 



Gene expression under water stress 

We initially compared each genotype under stress and 
control conditions to identify both general and genotype- 
specific stress-modulated genes, which were defined as 
upregulated (induced by water stress) or downregulated 
(repressed by water stress), as mentioned above. We iden- 
tified a total of 1817 probe sets that were modulated by 
water stress in at least one of the genotypes (Figure 3), in- 
cluding 1282 genes in the wild-type plants, 1422 genes in 
the mutant plants and 887 in both. 

We also compared the wild-type and mutant plants 
under water stress conditions, to identify genes that 
were differently-expressed under stress. We identified 
166 probe sets that were defined either as induced or re- 
pressed in the mutant plants, according to the level of 
expression under water stress conditions. 

Genotype-independent stress-response genes 

From the 887 probe sets modulated in both genotypes, 
286 were upregulated and 601 were downregulated 
under water stress (Figure 3). The upregulated genes pre- 
dominantly represented the following functional categor- 
ies: stress response, carbon compound and carbohydrate 



Table 3 Pigment content in leaves 



Pigment 






Control 








Water stress 




WT 


M 


IVl vs WT (%) 


P-value 


WT 


iVI 


M vs WT (%) 


P-value 


Chlorophyll a 


16308 


9704 


-40 


0.0315 


7607 


6804 


-11 


n.s. 


Chlorophyll b 


4700 


2374 


-49 


0.0106 


1901 


1609 


-15 


n.s. 


Total chlorophylls 


21007 


12078 


-43 


0.0247 


9508 


8413 


-12 


n.s. 


Ratio a/b 


346 


4.08 


18 


0.0001 


4.03 


4.24 


5 


n.s. 


Violaxanthin (V) 


968 


614 


-37 


0.0293 


535 


453 


-15 


n.s. 


Neoxanthin 


296 


171 


-42 


0.0095 


136 


97 


-28 


n.s. 


Antheraxanthin (A) 


147 


90 


-38 


n.s. 


86 


39 


-54 


0.0152 


Lutein 


1745 


856 


-51 


0.0071 


834 


649 


-22 


n.s. 


Zeaxanthin (Z) 


253 


258 


2 


n.s. 


224 


214 


-5 


n.s. 


P-Carotene 


1205 


734 


-39 


0.0464 


598 


566 


-5 


n.s. 


Total carotenoids 


4614 


2723 


-41 


0.0162 


2357 


2019 


-14 


n.s. 


(0.5A4-Z)/(VAZ) 


0.24 


0.32 


33 


0.0096 


0.34 


0.33 


-3 


n.s. 



Pigment quantification (|jg/g dw} of the flag leaf of wild-type (WT) and pale-green mutant (M) plants at anthesis under control and water stress conditions was 
performed by HPLC. Values are means of four samples. Xanthophylls de-epoxidation rate [(0.5A + Z)/(VAZ)]. 
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Figure 2 Principal component analysis plot of the array hybridization data. The x and y axis represent the two principal components of the 
total variance 59.31% and 15.76% respectively. The squares indicate the plants grown under control condition (CC), and the triangles indicate the 
plants grown under water stress (ES). The three biological replicates of each genotype: the pale-green mutant (IVl) and the wild-type (WT) are 
shown in the graph. 



Table 4 Significant functional categories under control conditions 



Regulation 


Functional categories 


Number of genes 


P-value 


Induced 


Ol.OSC-compound and carbohydrate metabolism 


4 


3.81 E-03 




01.05.03 polysaccharide metabolism 


3 


1.21 E-04 




01.05.03.03 chitin metabolism 


2 


1 .48E-05 




01.05.03.03.07 chitin catabolism 


2 


1.48E-05 




32 CELL RESCUE, DEFENSE AND VIRULENCE 


4 


2.25E-03 




36 SYSTEMIC INTERACTION WITH THE ENVIRONIVIENT 


3 


3.46E-03 




36.20 plant/fungal specific systemic sensing and response 


3 


2.69E-03 




36.20.16 plant defense response 


3 


1.73E-05 




36.20.16.03 jasmonic acid/ethylene dependent systemic resistance 


2 


488E-05 




36.20.16.05 systemic acquired resistance 


2 


9.10E-05 


Repressed 


01.03.16.03 DNA degradation 


1 


3.37E-03 




32. CELL RESCUE, DEFENSE AND VIRULENCE 


4 


2.25E-03 




32.01 stress response 


4 


2.99E-04 




34. INTERACTION WITH THE ENVIRONMENT 


4 


3.86E-03 




34.1 1 cellular sensing and response to external stimulus 


4 


2.65E-03 




34.11.03 chemoperception and response 


4 


3.12E-04 




34.1 1 .03.1 2 water response 


3 


1.53E-05 




36. SYSTEIVIIC INTERACTION WITH THE ENVIRONMENT 


5 


8.98E-06 




36.20 plant/fungal specific systemic sensing and response 


4 


1 .48E-04 




36.20.18 plant hormonal regulation 


4 


7.93 E-05 




36.20.18.05 abscisic acid response 


3 


3.78E-05 



Significant functional categories represented in the induced and repressed probe sets in the mutant compared to wild-type plants at anthesis under control condi- 
tions according to the MIPS Functional Catalogue Database [122]. Only functional categories with a cut-off of P< 0.005 were considered. 
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Table 5 Modulated probe sets in the mutant compared to the wild-type at anthesis under control conditions 



Expression 


Probe set 


Gene annotation 


Species^ 


UnlProt 


E- 
value 


Fold 
change 

(Wl vs WT) 


qRT-PCR 
fold 
expression 


nduced 


TaAffx.l28418.43.S1_at 


Chitinase 3 


; 


Q8W427 


1E-46 


5.8 


5 




Ta.28.1.S1_at 


Beta-l,3-glucanase 


; 


Q4JH28 


1E-133 


4.9 


- 




TaAffx.26815.1.S1_at 


Blufensini 


2 


B8X450 


1E-13 


44 


34 




Ta.169.1.S1_x_at 


Germin-like protein (GLP2b) 


1 


Q9SM34 


1E-154 


4.3 


2.8 




Ta.20188.1.S1_a_at 


Unknown 








4 


_ 




TaAffx.6092.1.Sl_at 


Naringenin,2-oxoglutarate 3-dioxygenase 


3 


B6SZR4 


1E-126 


3.6 


3.5 




Ta.3647.1.S1_at 


Putative uncliaracterized protein Sb01g018690 


4 


G5WXG4 


6E-51 


3.2 


- 




TaAffx.987441.S1_at 


Putative uncliaracterized protein Sb05g022350 


4 


C5Y4X4 


3E-60 


3.2 


- 




Ta.15159.1.S1_at 


Pathogenesis-related protein precursor 


2 


P93181 


3E-05 


3.1 


- 




Ta.20188.1.S1_x_at 


Unknown 








3.1 


- 




TaAffx.24475.1.S1_x_at 


Glucan endo-1,3-beta-glucosidase Gii 


2 


P15737 


9E-19 


3.1 


- 




Ta.22572.1.S1_at 


Unknown 








2.7 


- 




Ta.8614.2.S1_x_at 


WRKY45-like transcription factor 


1 


F8WPI8 


2E-53 


2.7 


2.2 




Ta.3830.1.Sl_a_at 


Unknown 








2.6 


- 




Ta.30028.1.Sl_s_at 


Peroxidase 10 (POXIO) 


5 


Q5I3E8 


1E-66 


24 






Ta.3830.2.Sl_x_at 


Unknown 








23 


- 




TaAffx.119315.2.S1_x_at 


Beta-glucanase 


2 


07M1K2 


5E-72 


2.2 


- 




Ta.2278.3.S1_x_at 


Ghitinase li (GHT2) 


1 


Q9XEN3 


3E-74 


2.2 


- 




Ta.3830.3.Sl_x_at 


Unknown 








2.1 


- 




TaAffx.24475.1.S1_at 


Glucan endo-1,3-beta-glucosidase Gil 


2 


PI 5737 


9E-19 


21 


- 




Ta.l6980.1.Sl_at 


Peroxidase 10 (POXIO) 


5 


Q5I3E8 


5E-10 


21 


- 




Ta.l6599.1.S1_x_at 


Adhesive/proline-rich protein 


3 


B4FH66 


8E-29 


2.1 


- 




Ta.8228.1.S1_at 


Agmatine coumaroyi transferase-1 (ACT-1) 


2 


A9ZPJ6 


0 


2 


2.5 


Repressed 


Ta.2638.1.S1_at 


RAB protein 


; 


041579 


1E-175 


-3.9 


0.08 




Ta.13255.1.S1_at 


Dehydrin (WZYl-1) 




Q8W192 


6E-93 


-3.9 


0.3 




Ta.23797.1.S1_x_at 


Late embryogenesis abundant protein (LEA1) 


, 


Q8GV49 


1E-104 


-3.6 


- 




Ta.24453.1.S1_5_at 


Oxalate oxidase (Germin protein) 


; 


09LD27 


7E-39 


-3.3 


- 




TaAffx.983941.S1_at 


0-type C2H2 zinc finger protein (ZFP23) 


; 


B0ZYY9 


1E-97 


-3.3 


0.7 




Ta.10390.1.S1_at 


Papain-like cysteine proteinase (PAP-14) 


2 


B4ESF2 


0 


-2.9 


- 




Ta.5497.1.Al_at 


Endonuclease 


2 


081958 


3E-13 


-2.8 


- 




Ta.6015.1.S1_at 


TdS40 protein (Fragment), putative 


6 


C4WYH7 


5E-83 


-2.7 


- 




TaAffx.83636.1.S1_at 


Unknown 








-2.7 






Ta4n0.1.S1_at 


lndole-3-acetic acid-amido synthetase (GH3.2), probable 


7 


P0G0M2 


1E-124 


-24 






Ta.5879.1.Sl_at 


Peptide chain release factor subunit 1 (ERFl), putative 


7 


Q75K79 


2E-92 


-24 






Ta.l 3682.1. Al_at 


Progesterone 5-beta-reductase, putative 


7 


O60DC5 


1E-35 


-2.3 






Ta.7543.1.S1_at 


Unknown 








-2.2 






Ta.5497.1.A1_x_at 


Endonuclease 


2 


081958 


3E-13 


-2.1 






Ta.l2477.1.S1_at 


Homeobox-leucine zipper protein (ATHB-6) 


3 


B6U539 


3E-03 


-2.1 





^Species: 1: Triticum aestivum, 2: Hordeum vulgare, 3: Zea mays, 4: Sorghum vulgare, 5: Triiicum monococcum, 6: Triticum turgidum ssp. Durum, 7: Oryza saiiva. 
-: not determined. 



metabolism, metabolite transport, amino acid metabol- with water stress are shown in the Additional file 3: 
ism, ABA response and proline biosynthesis (Additional Table S2. These include genes encoding sucrose synthase 
file 2: Table SI). Genes already known to be associated (SUSS, SUS4), glucose-6-phosphate dehydrogenase (G6PD4), 
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Up-regulated probe sets 



Down-regulated probe sets 



M Co-regulated 
157 286 





1 


M 




378 




Co-regulated 
601 



WT 
13d 



Figure 3 Modulated probe-sets. Venn diagrams show tlie number of modulated probe sets by the wild-type (WT) and mutant (IVl) at anthesis 
under water stress conditions compared with the control. 



branched-chain-amino-acid aminotransferase (BCAT3) 
and branched-chain alpha-keto acid dehydrogenase El 
a subunit (BCKDHEla), delta l-pyrroline-5-carboxyl- 
ate synthetase 2 (P5CS2), pyrroline-5-carboxylate re- 
ductase (P5CR) and RAB18. The downregulated genes 
predominantly represented the following functional 
categories: phosphate and secondary metabolism, e.g. 1- 
aminocyclopropane-l-carboxylate synthase (ACS) and 12- 
oxophytodienoate reductase (OPR); interaction with the 
environment, e.g. chitinases, |3-glucanases and pathogenesis- 
related proteins; immune response, e.g. transcription factor 
WRKY; and phenylpropanoid metabolism, e.g. phenylal- 
anine ammonia-lyase. These genotype-independent stress- 
modulated genes correspond to the well-characterized 
drought stress response of wheat, confirming that the mu- 
tant plants preserve much of the stress-response signaling 
network that allows wheat plants to tolerate water stress 
[34,43]. 

Cenotype-dependent stress-response genes in wild-type 
plants 

We identified 395 probe sets that were modulated solely 
in the wild-type plants or to a greater extent (above the 2- 
fold threshold) than in the mutant plants. This included 
265 upregulated and 130 downregulated genes (Figure 3). 
The upregulated genes predominantly represented carbo- 
hydrate and pyruvate metabolism and the glyoxylate cycle, 
whereas the downregulated represented metabolism, inter- 
action with the environment, protein modification and sig- 
naling pathways (Table 6). Selected genes specifically 
modulated by the wild-type plants under stress are listed 
in Table 7, including carbon fixation enzymes such as mal- 
ate dehydrogenase (NADP-ME3) and phosphoenolpyr- 
uvate caboxykinase (PCKl), and glyoxylate cycle enzymes 
such as isocitrate lyase (ICL) and malate synthase (MLS). 

Genotype-dependent stress-response genes in 
mutant plants 

We identified 535 probe sets that were modulated solely in 
the mutant plants or to a greater extent (above the 2-fold 



threshold) than in the wild- type plants. This included 157 
upregulated and 378 downregulated genes representing di- 
verse functional categories (Figure 3; Table 8). The most 
significant categories in the upregulated group were associ- 
ated with chloroplast transcription, interactions with the 
environment (cellular sensing, response to external stimuli, 
chemoperception and response), stress responses and me- 
tabolism, whereas the downregulated group was dominated 
by functional categories such as interaction with the envir- 
onment (response to external stimuli), stress responses 
and plant defense (Table 8). The downregulated genes also 
included those induced in the mutant under control 
conditions. 

The genes that were upregulated specifically in the 
mutant included those involved in the synthesis of caroten- 
oid precursors, e.g. hydroxymethylglutaryl-CoA synthase 
(MVAl), l-deoxy-D-xylulose-5-phosphate synthase (DXS) 
and phytoene synthase (PSY), and those involved in 
chloroplast gene expression, photosynthesis, chlorophyll 
binding and detoxification, e.g. glutamyl-tRNA reductase 1 
(HEMAl), light harvesting complex II components such as 
LHB1B2 and ATPase subunit delta, and the Rubisco small 
subunit. We observed the downregulation of several genes 
representing glutamate metabolism such as glutamate de- 
carboxylase (GAD), glutamate dehydrogenase (GDH), y- 
glutamyltranspeptidase (GGT) and L-aspartate oxidase 
(AO); and several signaling pathway components, including 
calcineurin B-like protein 1 (CBLl), inositol phosphatase 
(SALl) and MPK3 (Table 7). 

Differentially-expressed stress-response genes under 
stress conditions 

We compared the transcriptional profiles of wOd-type and 
mutant plants under water stress and identified 166 
differently-expressed probe sets (29 induced and 137 re- 
pressed in the mutant). The only functional category signifi- 
cantly represented among the 29 induced probe sets 
was the channel/pore transport class (Additional file 4: 
Table S3), corresponding to the transcripts for two plasma 
membrane intrinsic proteins (PIP2;2, PIP2;5) that were 
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Table 6 Significant functional categories in the wild-type under water stress 



Regulation 



Functional categories 


Number of genes 


P-value 


01 METABOLISM 


53 


6.57E-11 


01.01.11 metabolism of the pyruvate family (alanine, isoleucine, leucine, valine) and D-alanine 


3 


1 .90E-03 


01.01.1 1.04 metabolism of leucine 


3 


3.37E-04 


01. 01. 11. 04.02 degradation of leucine 


3 


1.49E-05 


01.03.16.03 DNA degradation 


2 


8.84E-04 


Ol.OSC-compound and carbohydrate metabolism 


33 


211E-10 


01.05.02 sugar, glucoside, polyol and carboxylate metabolism 


21 


1 .96E-08 


01.05.03 polysaccharide metabolism 


9 


9.71 E-06 


01.05.05C-1 compound metabolism 


5 


435E-03 


01.05.07.04C-3 compound anabolism 


2 


3.23E-03 


02 ENERGY 


8 


4.60E-03 


02.04 glyoxylate cycle 


4 


2.08E-07 


36.25 animal specific systemic sensing and response 


4 


3.33E-03 


36.25.16 immune response 


4 


2.61 E-03 


42.01 cell wall 


6 


4.09E-03 


70.10.03 chromosome 


7 


1.71 E-05 


01 METABOLISM 


33 


3.37E-09 


01.04 phosphate metabolism 


17 


1.29E-06 


01.05C-compound and carbohydrate metabolism 


11 


1.73E-03 


01.20 secondary metabolism 


5 


3.56E-03 


01.20.35 metabolism of secondary products derived from L-phenylalanine and L-tyrosine 


3 


3.70E-03 


01.20.35.01 metabolism of phenylpropanoids 


3 


2.81 E-03 


14 PROTEIN FATE (folding, modification, destination) 


16 


2.28E-03 


14.07 protein modification 


14 


5.23E-04 


14.07.03 modification by phosphorylation, dephosphorylation, autophosphorylation 


12 


3.51 E-05 


20.01.01.07 anion transport 


3 


1 .90E-03 


30.01.05 enzyme mediated signal transduction 


6 


6.96E-04 


30.01.05.01 protein kinase 


6 


5.69E-05 


30.01.05.01.06 serine/threonine kinase 


5 


8.08E-05 


30.05 transmembrane signal transduction 


5 


1.14E-03 


30.05.01 receptor enzyme mediated signalling 


5 


5.10E-04 


30.05.01.12 transmembrane receptor protein tyrosine kinase signalling pathways 


5 


1.54E-05 


32 CELL RESCUE, DEFENSE AND VIRULENCE 


10 


2.00E-03 


32.01 stress response 


9 


1 .59E-04 


32.01.06 cold shock response 


4 


6.69E-04 


34 INTERACTION WITH THE ENVIRONMENT 


11 


1.77E-03 


34.1 1 cellular sensing and response to external stimulus 


11 


7.70E-04 


34.1 1 .03.1 2 water response 


3 


2.94E-03 


34.1 1.09 temperature perception and response 


5 


4.24E-04 


70.03 cytoplasm 


8 


3.86E-03 



Upregulated 



Downregulated 



Significant functional categories represented in the specifically upregulated and downregulated probe sets in wild-type plants at anthesis under water stress 
conditions compared to well-watered (control) conditions according to the MIPS Functional Catalogue Database [122]. Only functional categories with a cut-off of 
P< 0.005 were considered. 
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Table 7 Selected genotype-dependent regulated genes under water stress 



Probe set ID 



Gene name 



Arabidopsis ID Log2 ratio WT Log2 ratio M Log2 ratio M/WT 

stress/WT control stress/M control under water stres 



PHOTOSYNTHESIS AND PIGMENT SYNTHESIS 

Ta.23025.1.A1_at 
Ta.28496.1.A1_x_at 



LHB1B2; Light-harvesting 
complex II chlorophyll 
a/b binding protein 1 



Ta.30702.1.Sl_x_at 

Ta.3249.2.Sl_x_at 

Ta.3249.3.A1_at 

Ta.l 8362.1. A1_at 

Ta.23419.1.Sl_x_at 

Ta.8718.2.Sl_at 

Ta.8718.2.Sl_x_at 

Ta.3243.1.Sl_at 

Ta.3243.1.Sl_x_at 

Ta.20830.1.A1_at 

TaAffx.128414.24.A1_s_at 
Ta.l 7243.1. SI at 



Ta.3366.1.S1_at 



Ta.20776.1.S1 at 



ELIPl; Early light-inducible 
protein 

ATPD; F-type H* -transporting 
ATPase subunit delta 

HEIVlAl; Glutamyl- 
tRNAreductase 1 

RUBISCO small subunit; Ribulose 
bisphosphate carboxylase small 
chains, chloroplast precursor 

RUBISCO large subunit 

MVAl; Acetyl-CoA C- 

acetyltransferase/ 

hyd roxymethylg I utary l-CoA 

synthase 

DXS; 1-deoxy-D-xylulose-5- 
phosphate synthase; CLAl 
(Cloroplastos Alterados 1); 

PSY; Phytoene synthase 



OXIDATIVE STRESS AND ROS SIGNALING 
Detoxification 



Ta.11386.2.Sl_a_at 

Ta.l1386.2.Sl_x_at 

Ta.488.2.S1_at 

Ta.23079.1.Sl_a_at 

Ta.30028.1.Sl_s_at 

TaAffx.613.2.A1_x_at 

Ta.l 8560.1. Sl_at 

Ta.21 505.1. Sl_at 

Ta.30697.1.Sl_at 

TaAffx.39568.2.Sl_at 

Ta.6127.1.Sl_at 

TaAffx.l 00459. 1.Sl_at 

Ta.8571.1.S1_x_at 

TaAffx.79142.1.S1_at 

Ta.8571.1.Sl_a_at 
Ta.8571.1.S1_x_at 

Ta.237041.S1_s_at 
Ta.303.2.Sl_at 



TAPX; Thylakoidal ascorbate 
peroxidase 

APX4; Ascorbate peroxidase 4 
APXl; Ascorbate peroxidase 1 
PERI 2; Peroxidase 12 
Peroxidase, putative 



Peroxiredoxin type 2. putative 

Glutaredoxin family protein 

Glyoxalase I family protein 

ATGLXl; Glyoxalase I homolog; 
lactoylglutathione lyase 

Glyoxalase I family protein; 
lactoylglutathione lyase 
family protein 

ATGSTF13; Glutathione 
transferase 



AT2G34420 



AT3G22840 
AT4G09650 
ATI G58290 
TC371960'' 

ATCG00490 
AT4G 11820 

AT4G 15550 
AT5G 17230 

ATI G77490 

AT4G09010 
ATI G07890 
AT1G71695 
AT4G31760 
AT4G33420 
ATI G49570 
AT5G05340 
AT5G58390 
AT3G52960 
AT3G62950 
AT1G80160 
ATI Gl 1840 

AT1G80160 
AT3G62760 



1.14 



-1.80 
-1.37 



1.12 



1.22 
1.32 



2.04 
1.63 
1.13 
1.20 



1.47 
1.29 
1.52 
1.47 
1.36 



1.13 



1.55 



1.07 



1.25 
1.31 
1.7 
-1.03 
-1.62 

-1.78 
-1.29 
-1.05 
-1.57 



-1.03 
-0.76 

-1.04 
-1.03 

-1.36 
-1.08 



1.02 



1.41 
0.95 



1.04 



-0.85 



1.73 



1.07 



-1.07 
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Table 7 Selected genotype-dependent regulated genes under water stress (Continued) 



Ta.3681.1.Sl_x_at 
TaAffx.79276.1.Sl_x_at 

Ta.21001.1.S1_at 
TaAffx.64682.1.Sl_at 

Signaling 

Ta.4936.1.Sl_at 

Ta.8465.1.Sl_at 

Ta.5041.1.Sl_at 
Ta.5041.2.Sl_a_at 
Ta.5272.3.Sl_x_at 
Ta.236.1.Sl at 



ATGSTU 19; Glutathione 
S-transferase TAU 19 

ERD9; Early-responsive to 
dehydration 9; glutathione 
transferase 

GGTl; Gamma-glutamyl 
transpeptidase 1 

RBOHD; Respiratory burst 
oxidase homologue D 



CIPK2; CBL-lnteracting protein 
kinase 2 

CIPK5; CBL-lnteracting protein 
kinase 6 

CIPK23; CBL-lnteracting protein 
kinase 23 

CBLl; Calcineurin B-like protein 
MAPK3; MAP kinase 3 



CARBOHYDRATE METABOLISM 

Galactose and Raffinose family metabolism 



Ta.26092.1.Sl_at 

Ta.6023.1.Sl_at 

Ta.75412.Sl_a_at 
Trehalose synthesis 

Ta.20649.1.Sl_x_at 

Starch degradation 

Ta.10281.1.Sl_a_at 

Ta.l6135.1.A1_at 

TaAffx.68872.1.S1_at 

Ta.8827.1.Sl_at 

Sucrose hydrolisis 

Ta.9000.1.Sl_at 

Ta.2788.1.Al_at 

Ta.27892.Sl_at 

Cellulose synthesis 

Ta.1 2382.1. Sl_at 
Ta.1 2382.1. Sl_x_at 
Ta.7192.1.S1_at 

Ta.4447.1.S1_at 
Ta.4447.1.S1_x_at 
Ta.4447.2.S1 x at 



SIP1; Putative galactinol-sucrose 
galactosyl transferase 1; 
Raffinose synthase 

SIP2; Putative galactinol-sucrose 
galactosyl transferase 2 

GalAK; Galactokinase 



ATTPSe; Alpha.alpha-trehalose- 
phosphate synthase 6 



AMY1; Alpha-amylase-like 

CT-BMY; Chloroplast 
Beta-Amylase 



AT1G78380 
ATI G 10370 

AT4G39540 
AT5G47910 

AT5G07070 
AT4G30960 
AT1G30270 

AT4G17615 
AT3G45640 

AT1G55740 

AT3G57520 
AT3G 10700 

ATI G68020 



AT4G25000 
AT4G 17090 



-1.02 



DPE2; Disproportionating enzyme 2 AT2G40840 



BFRUCT3; Beta-fructosidase; ATI G52660 

invertase, vacuolar 

Sucrose:sucrose 1- B5TK35'' 
fructosyltransferase 

Sucrose:fructan 5- 096466" 
fructosyltransferase 



ATCSLB02; Cellulose synthase AT2G32620 

ATCSLCl 2; Cellulose synthase AT4G07960 
like CI 2 

CESAl; Cellulose synthase 1 AT4G32410 

CESA3; Cellulose synthase 3 AT5G05170 



1.23 



1.46 



1.06 



1.20 



1.27 



1.47 



1.21 

-1.51 
-1.25 



1.51 



1.04 
1.11 
1.15 

1.24 
1.40 
1.16 



-1.23 

-0.99 
1.07 



.59 
.22 
.04 
.20 



1.06 



1.14 



1.09 



1.48 



1.49 



1.33 
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Table 7 Selected genotype-dependent regulated genes under water stress (Continued) 



Ta.l2382.3.Sl_at 
Ta.5628.1.A1_at 



ATCSLB05; Cellulose synthase 

IRX9; Irregular xylem 9; 
xylosyltransferase 

XT2; UDP-Xylosyltransferase 2 

Xyloglucan:xyloglucosyl 
transferase, putative 



Ta.29111.1.A1_at 
Ta.l 3337.1. Sl_at 
Ta.13337.2.Sl_at 
Ta.13337.2.Sl_x_at 
TCA AND RELATED PATHWAYS 
Pyruvate metabolism and glyoxylate cycle 



Ta.l 8775. l.Sl_at 
Ta.25543.2.Sl_at 

Ta.l 3280.1. Sl_a_at 

Ta.l3280.2.Al_at 

Ta.23970.1.Al_x_at 

Ta.23989.1.A1_at 

TaAffx.79218.1.Sl_at 

Ta.l 9423.1. Sl_a_at 

Ta.l 3281. l.Sl_at 

Amino acid metabolism 

Ta.l5150.1.Sl_at 

Ta.l5150.1.Sl_x_at 

Ta.30569.1.Sl_a_at 

Ta.l 2252.1. Sl_s_at 

Glutamate metabolism 

TaAffx.l 29066.1. Sl_at 
Ta.25990.1.Al_at 
Ta.25990.1.Al_x_at 
Ta.l 870.1. Sl_a_at 
Ta.21001.1.Sl_at 

Ta.30684.1.Sl at 



NADP-ME4; NADP-malic 
enzyme 4 

NADP-ME3; Malate 
dehydrogenase (oxaloacetate- 
decarboxylating) 

PCKl; Phosphoenolpyruvate 
carboxykinase 

MLS; Malate synthase 
ICL; Isocitrate lyase 

PDC2; Pyruvate decarboxylase-2 
Pyruvate decarboxylase, putative 



MCCA; Methylcrotonoyl-CoA 
carboxylase subunit alpha 

MCCB; Methylcrotonoyl-CoA 
carboxylase beta chain 

BCE2; 2-oxoisovalerate 
dehydrogenase E2 component 
(dihydrolipoyltransacylase) 



GAD; Glutamate decarboxylase 1 
GAD3; Glutamate decarboxylase 3 

GDH2; Glutamate dehydrogenase 2 

GGTl; Gamma-glutamyl 
transpeptidase 1 

AO; L-aspartate oxidase 



REGULATORS OF CARBON METABOLISIVl AND SENSING 



Ta.l 2733.1. Sl_at 

Ta.4492.1.Sl_at 

Ta.l 01 30.1. Sl_at 
Ta.l 3336.1. Sl_at 
Ta.21 035.1. Sl_at 
Ta.27316.1.Sl_at 
TaAlfx.l8447.1.Sl_at 
TaAf[>c.l8447.3.Sl_s_at 
TaAffx.l84475.Sl s at 



SALl; Inositol or 

phosphatidylinositol phosphatase 

PI3P5K; Phosphatidylinositol-3- 
phosphate 5-kinase 

RAP2.12; Related to AP2 12 

RAP2.7; Related to AP2 7 

RAP2.8; RAV2 (Regulator of the 
ATPase of the vacuolar membrane) 

BZIP63; Basic Leucine 
Zipper 63; BZ02H3 



AT4G 15290 
AT2G37090 

AT4G02500 
AT5G57540 



1.09 



AT1G79750 
AT5G25880 

AT4G37870 

AT5G03860 
AT3G21720 

AT5G54960 
AT5G01320 

ATI G03090 

AT4G34030 
AT3G06850 



AT5G 17330 
AT2G02000 

AT5G07440 
AT4G39640 

AT5G 14760 

AT5G63980 

AT1G71010 

AT1G53910 
AT2G28550 
ATI G68840 

AT5G28770 



1.81 

1.26 
1.38 
1.21 
1.44 



1.11 
1.20 

1.05 
1.10 
1.92 
2.38 
1.63 



1.09 
1.12 
1.10 

1.13 



1.05 
1.68 



1.14 
1.07 



-1.42 
-1.53 
-1.31 
-0.99 



-1.02 
1.02 



-2.02 
-1.11 
-1.84 
-2.16 
-1.55 



-1.30 

-1.32 
-1.59 
-1.41 
-1.57 



-1.31 
-1.90 
-1.35 



-1.05 
-1.11 
-1.23 



-1.32 
-1.05 
-1.44 

-1.05 
-1.42 
-1.03 
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Table 7 Selected genotype-dependent regulated genes under water stress (Continued) 



CHLOROPLAST TRANSCRIPTION AND TRANSLATION 



TaAffx.271S6.1.Sl_at 
TaAffx.1 07547. 1.Sl_at 
TaAffx.11 3923.1. Sl_at 

Ta.769.1.S1_at 

TaAffx.1 28946.4Sl_at 
OTHERS 

Ta.51613.Sl_at 
Ta.51613.Sl_x_at 
Ta.6870.2.Sl a at 



RNA polymerase beta' subunit-2 

RNA polymerase beta' subunit-1 

Chloroplast DNA-dependent RNA 
polymerase B subunit 

PTAC16; Plastid transcriptionally 
active 16 

Chloroplast ribosomal protein S2 



CRK29; Cystein-rich RLK 
(Receptor-like protein kinase) 29 

MKK6; MAP kinase kinase 6 



ATCG00170 
ATCG00180 
ATCG00190 

AT3G46780 

ATCG00160 

AT4G21410 

AT5G56580 



1.24 
1.08 

1.26 

1.08 
1.18 

-1.73 
-132 



Selected genes were specifically significant regulated by the wild-type (WT) and the mutant (M) under water stress compared with the same genotype under well- 
watered (control) conditions. Last column (M vs WT) compares mutant versus WT expression under water stress. Positive and negative values indicate upregulation 
and downregulation respectively, whereas no statistically significant differences are indicated with a dash (-). 
^UniProt ID from Triticum turgidum subsp. durum. 

''Blast executed against TaGI (DFCI) Wheat Gene Index: TC371960. homologue to UniRef100_P26667. 



downregulated in both genotypes but to a lesser extent 
in the mutant (Additional file 3: Table S2). Other genes 
that were induced in the mutant under water stress in- 
cluded those encoding ascorbate peroxidase (APX4), su- 
crose: fructan 6-fructosyltransferase (SFT6), mjo-inositol 
oxygenase (MIOXl) and dihydrodipicolinate reductase 3 
(crrl) (Table 7 and Additional file 3: Table S2). 

The genes that were repressed in the mutant plants under 
stress were predominantiy related to carbohydrate metabo- 
lism, peroxidase/oxidative stress responses, the glyoxylate 
cycle and chromosome subcellular localization (Additional 
file 4: Table S3). Representative repressed genes included 
those encoding RAP2.8, GDH, GAD, CBL-interacting pro- 
tein kinase 23 (CIPK23) and glyoxalase (Table 7). We also 
found genes encoding ICL, MLS and trehalose phosphate- 
synthase (TPS) which were not affected by stress in mutant 
plants but were upregulated in wild-type plants, thus 
appearing in the repressed transcript category (Table 7). 

Proteomic profiling 

The effect of the mutation was also investigated by car- 
rying out a comparative proteomic analysis between 
wild-type and mutant plants at anthesis under the same 
conditions as described above. 

When wild-type and mutant plants were compared under 
well-watered (control) conditions, we matched 182 protein 
spots (data not shown) among which two showed a signifi- 
cant (>2-fold) difference in signal intensity between geno- 
types, although neither spot could be identified because of 
the low abundance of the corresponding protein. Even so, 
this experiment supported the transcriptomic and sug- 
gested that there were only minor differences between the 
genotypes under normal conditions. 

When we compared each genotype under stress and 
control conditions, we matched 207 protein spots 



among both genotypes. From the 207, we found that 32 
spots were modulated by water stress in the wild-type 
plants and 62 were modulated in the mutant plants, in- 
cluding 25 spots that were common to both genotypes. 
This left seven modulated spots specific to the wild-type 
plants and 37 specific to the mutant. 

When the two genotypes were compared under water 
stress, we identified 28 spots that were more abundant 
in the mutant and 17 that were less abundant in the 
mutant. We selected 19 of the induced spots and all 17 
of the repressed spots based on the likelihood of extract- 
ing useful amounts of protein, and we identified 13 pro- 
teins from the induced set and 8 from the repressed set 
(Additional file 5: Figure S2). The proteins were identi- 
fied by peptide mass fingerprinting (Additional file 6: 
Table S4) and MS/MS analysis (Additional file 7: Table S5). 
The induced proteins predominantiy represented carbohy- 
drate and amino acid metabolism, photosynthesis, de- 
toxification and phytoalexin biosynthesis, whereas the 
repressed proteins predominantly represented photo- 
synthesis (Table 9). Two proteins that were specifically 
upregulated in the mutant were the oxygen-evolving 
enhancer protein 2 corresponding to the PSBP subunit 
of photosystem II (A.n. Q00434), and catalase-1 (CATAl; 
A.n. Q00434). Two proteins that were specifically down- 
regulated in the mutant were the Rubisco small subunit 
chain PW9 (A.n. P26667) and the photosystem I reaction 
center subunit IV (A.n. P13194). 

Discussion 

Chloroplast ultrastructure and pigment content 

Chlorotic mutants with affected chloroplast ultrastructure 
have been described in several reports [9,14,44-46]. The 
chloroplasts in these plants often contain fewer grana 
structures than wild-type plants, and thinner thylakoid 
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Table 8 Significant functional categories in the mutant under water stress 



Regulation 



Functional categories 


Number of genes 


P-value 


01 METABOLISM 


26 


3.94E-03 


11.02.03.01 general transcription activities 


3 


8.94E-04 


1 1.02.03.01.04 transcription elongation 


3 


3.51 E-06 


32 CELL RESCUE. DEFENSE AND VIRULENCE 


12 


1.61 E-03 


32.01 stress response 


11 


5.80E-05 


32.01.03 osmotic and salt stress response 


4 


3.98E-03 


34 INTERACTION WITH THE ENVIRONMENT 


16 


4.91 E-05 


34.1 1 cellular sensing and response to external stimulus 


14 


2.11E-04 


34.11.03 chemoperception and response 


8 


4.89E-03 


34.11.03.13 osmosensing and response 


4 


419E-03 


36.20.18 plant hormonal regulation 


7 


2.53E-03 


36.20.18.02 ethylen response 


3 


3.86E-03 


01 METABOLISM 


72 


1.44E-09 


01.04 phosphate metabolism 


28 


2.37E-04 


Ol.OSC-compound and carbohydrate metabolism 


24 


1.03 E-03 


01.05.03.03 chitin metabolism 


3 


1.33E-04 


01.05.03.03.07 chitin catabolism 


3 


1.33E-04 


01.20 secondary metabolism 


11 


3.27E-04 


01.20.35 metabolism of secondary products derived from L-phenylalanine and L-tyrosine 


5 


2.65E-03 


01.20.35.01 metabolism of phenylpropanoids 


5 


1.73E-03 


01.20.35.01.03 metabolism of lignins 


4 


3.58E-04 


14.07.03 modification by phosphorylation, dephosphorylation. autophosphorylation 


20 


8.53E-04 


16.13C-compound binding 


4 


2.49E-03 


16.17.01 calcium binding 


7 


2.97E-04 


16.25 oxygen binding 


8 


3.57E-04 


20 CELLULAR TRANSPORT TRANSPORT FACILITIES AND TRANSPORT ROUTES 


37 


1 .41 E-05 


20.01 transported compounds (substrates) 


29 


1 .30E-04 


32 CELL RESCUE. DEFENSE AND VIRULENCE 


29 


6.36E-07 


32.01 stress response 


20 


3.59E-06 


32.01.01 oxidative stress response 


9 


1.79E-05 


32.01.06 cold shock response 


7 


2.23 E-04 


32.07 detoxification 


8 


8.06E-04 


32.07.07 oxygen and radical detoxification 


8 


7.1 3 E-04 


32.07.07.05 peroxidase reaction 


5 


6.49E-04 


34 INTERACTION WITH THE ENVIRONMENT 


35 


1.35E-08 


34.1 1 cellular sensing and response to external stimulus 


34 


3.56E-09 


34.11.03 chemoperception and response 


16 


5.00E-04 


34.11.05 mechanical stimulus perception and response 


2 


1 .86E-03 


34.1 1.09 temperature perception and response 


10 


2.62E-05 


34.11.10 response to biotic stimulus 


14 


446E-08 


36 SYSTEMIC INTERACTION WITH THE ENVIRONMENT 


19 


3.43 E-06 


36.20 plant/fungal specific systemic sensing and response 


18 


3.84E-06 


36.20.16 plant defense response 


7 


3.72E-05 


36.20.16.05 systemic acquired resistance 


3 


1 .96E-03 



Upregulated 



Downregulated 
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Table 8 Significant functional categories in the mutant under water stress (Continued) 



36.20.18 plant hormonal regulation 


12 


1 .42E-03 


40.10 cell death 


6 


1.02E-04 


40.10.02 apoptosis (type 1 programmed cell death) 


6 


9.85 E-06 


40.10.02.01 anti-apoptosis 


6 


7.34E-08 


70.02 eukaryotic plasma membrane/membrane attached 


12 


4.61 E-04 



Significant functional categories represented in the specifically upregulated and downregulated probe sets lists in the mutant at anthesis under water stress 
conditions compared to well-watered (control) conditions according to the MIPS Functional Catalogue Database [1 22]. Only functional categories with a cut-off of 
P< 0.005 were considered. 



membranes that are not stacked in the normal manner 
[47]. In contrast, the chloroplasts in our pale-green mutant 
contained abundant grana with normally-stacked thyla- 
koids, but their distribution within the stroma was disorga- 
nized, potentially reflecting a lack of membrane continuity 
or the degradation of inter-grana lamella membranes. The 
stacking adhesion of the thylakoids in grana structures is 
regulated by the phosphorylation of photosystem II pro- 
teins, which also ensures photosynthetic efficiency [45]. We 
found no significant difference in the photosynthetic rate of 
the wild-type and mutant and no quantitative differences in 
transcripts or proteins representing the photosynthetic ap- 
paratus at anthesis under control conditions. 

The mutant chloroplasts also contained more stromal 
plastoglobules than wild-type chloroplasts. These lipopro- 
tein bodies are structurally continuous with the thylakoid 
membranes and facilitate content exchange with the 
thylakoid lumen [48]. They may also be involved in 
thylakoid development and disassembly, increasing in 
abundance in response to oxidative stress and deliver- 
ing the antioxidant molecules stored on them to the 
thylakoid lumen, where they scavenge ROS [49,50] . The 
structure of the inter-grana lamella in the mutant chlo- 
roplasts could therefore reflect thylakoid reorganization 
or the impact of oxidative stress already at tillering 
stage under field capacity conditions. 

Although the pigment content of the mutant plants 
was -41% lower than in wild-type plants under control 
conditions, the photosynthetic rate was similar (Table 1). 
This seems illogical given that carotenoids and chloro- 
phylls are required for the proper folding, assembly and 
stability of LHC apoproteins [51]. However, we detected 
a greater xanthophyll de-epoxidation rate in the mutant 
as indicated by the relationship (0.5A -i- Z)/(VAZ). The 
zeaxanthin in the mutant may compensate for the lack 
of other xanthophylls in the antenna proteins [52], and 
may also preserve the structure of the photosystem and 
its membranes by dissipating photochemical energy as 
heat to prevent photo-oxidation [22]. This may explain 
the abnormal chloroplast ultrastructure in the mutant 
but the near-normal photosynthetic rate. Our data 
suggest that, under control conditions, the mutant com- 
pensates for the lower pigment content and that minor 



oxidative stress has no detrimental effect on photosyn- 
thetic capacity or biomass accumulation. 

However, under water stress conditions the level of de- 
epoxidation increases in wild-type plants, and both geno- 
types show similar values for the relationship (0.5A -i- Z)/ 
(VAZ) (Table 3). This suggests that the mutant has a lower 
capacity for increasing the de-epoxidation rate. An alterna- 
tive explanation is that the pool of convertible violaxanthin 
is rapidly depleted and only non-convertible violaxanthin 
remains to maintain the structural integrity of the antenna 
[53]. Because xanthophylls help to prevent photo-oxidation, 
the mutant plants without additional de-epoxidated carot- 
enoids to cope for the water stress damage may sense 
higher oxidative stress in the chloroplast under water stress. 
Thus, the pale-green mutant may be more susceptible since 
the ability of plants to withstand oxidative stress determines 
the level of overall stress tolerance [54]. 

Differences between the wild-type and mutant 
transcriptomes at anthesis (no stress) 

There were only minor differences in gene expression be- 
tween the wild-type and mutant plants under control con- 
ditions at anthesis, mirroring the similar net photosynthetic 
rates (Table 1). The few repressed genes in the mutant were 
predominantly related to ABA responses and the induced 
genes mostly represented pathogenesis-related (PR) pro- 
teins and peroxidases involved in plant defense (Table 5). 

Apart from their role in defense, peroxidases may be 
modulated by cellular redox disturbances [55]. In such 
cases, nuclear gene expression can be regulated by ROS- 
dependent signal transduction or low-molecular-weight 
antioxidants produced in the chloroplast. For example, 
the expression of PR proteins and peroxidases increases 
in ascorbate-deficient mutants [56]. Since no lesion or 
pathogen attack was detected in our investigation, it is 
most likely that low-level of oxidative stress triggers 
defense gene expression in the pale-green mutant. Pro- 
teins involved in plant defense are also sensitive to hor- 
mones such as salicylic acid (SA), jasmonate (JA) and 
ethylene, but we did not find any differently-expressed 
genes related to these pathways in the mutant. Plant 
defense signaling pathways are also activated through 
WRKY transcription factors that bind to pathogen- 



Table 9 Differently-expressed proteins In the mutant (M) compared to the wild-type (WT) under water stress 



Expression Function 


Protein description 


SwissProt accession 
number 


Species'" Ratio WT stress/WT 
control 


Ratio M stress/IVl 
control 


Ratio M/WT under 
water stres 


Method" 


Induced Carbohydrate metabolism 


G3PX; Glyceraldehyde-3-phosphate 
dehydrogenase, cytosolic 


P26517 


1 




304.7 


PMF 




G3PX; Glyceraldehyde-3-phosphate 
dehydrogenase, cytosolic 


P26517 


1 


- 


2154 


PMF 




Glyceraldehyde-3-phosphate 
dehydrogenase A. chloroplastic 


P09315 


2 


170.3 


87.8 


MS/MS 




p -1- L," 1 1 4- 1^1 m/^ni,^ 

Fructose-bisphosphatealdolaseF2CR16 


BAJ85287^ 


1 




3.7 


MS/MS 


Amino acid metabolism 


Cysteine synthase 


P38076 


3 


- 


252.7 


MS/MS 




Chloroplast aspartate aminotransferase 


ACG59771 = 


3 8.7e"^^ 




101.1 


MS/MS 




Putative 3-beta hydroxysteroid 
dehydrogenase/isomerase protein 


BAJ86056'' 


1 


1345 


75.8 


MS/MS 


Photosynthesis 


ATPB; ATP synthase subunit beta, chloroplastic 


P20858 


3 


- 


114.2 


PMF 




PSBO; Oxygen-evolving enhancer protein 
1. chloroplastic (PSIl) 


040459 


4 


- 


2.5 


PMF 




PSBP; Oxygen-evolving enhancer protein 
2. chloroplastic (PSIl) 


000434 


3 


2.1 


2.1 


PMF 


Biosynthesis phytoalexins 


soflavone reductase 


B5M699 


1 


- 


2.6 


PMF 


Translation 


Elongation factor Tu. chloroplastic-like. predicted 


XP_003575279" 


5 


2.2 


2.5 


MS/MS 


Detoxification 


CATAl; Catalase-1 


043206 


3 


3.4 


2.5 


PMF 


Repressed Photosynthesis 


Chloroplast light-harvesting chlorophyll a/b 
binding protein (Lhc II type! CAB) 


ADL41158" 


3 0.4 


0.3 


0.5 


MS/MS 




Photosystem 1 reaction center subunit IV. 
chloroplastic (PSI RCsublV) 


P13194 


1 


0.5 


0.4 


MS/MS 




Precursor of CP29. core chlorophyll a/b binding 
(CAB) protein (PSIl Lhc4) 


CAA44777" 


1 595.9 


2655 


0.4 


MS/MS 




RBS2; Ribulose bisphosphate carboxylase small 
chain PW9. chloroplastic 


P26667 


3 


0.5 


0.4 


PMF 




Chlorophyll a-b binding protein CAB 1B-21 
(PSI Lhc 1) 


ACO06083'' 


3 




0.3 


MS/MS 




Chlorophyll a-b binding protein lB-21. 
chloroplastic (CAB 1B-21) 


09SDM1 


1 


2.5e-" 


0.002 


MS/MS 


Heme binding 


HBL2; Non-symbiotic hemoglobin 2 


024521 


6 


05 


0.5 


PMF 




C71AJ; Cytochrome P450 71A19 


O9T0K0 


6 0.3 


04 


0.4 


PMF 



^Accession number from NCBI. 

'^Species: 1 : Hordeum vulgare; 2: Zea mays; 3:Triticum aestivum; 4: Nicotiana tabacum; 5: Brachypodium distachyon; 6: Arabidopsis thaliano. 
^Method: PMF: Peptide Mass Fingerprinting; MS/MS: MALDI TOF/TOF. 
Ratios are only indicated for differences statistically significant. 
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response and wound-inducible promoter elements (W- 
boxes and GCC-like elements). The TaWRKY4S gene 
was induced in the mutant, and may help to activate 
some of the PR genes we also identified. This transcrip- 
tion factor confers resistance to Fusarium head blight 
(FHB) in wheat [57] and increases the expression of PR- 
1 and PR-2 genes in rice to confer resistance to Pseudo- 
monas syringae, salinity and drought tolerance, as well 
as reducing ABA sensitivity [58]. 

Genes encoding defense enzymes that synthesize second- 
ary metabolites such as phenylpropanoids and flavonoids 
were also induced in the mutant, including naringerin, 2- 
oxoglutarate 3-dioxygenase and agmatine coumaroyltrans- 
ferase (ACT), the latter representing the final step in the 
hydroxycinnamic acid amide (HCAA) synthesis pathway. 
Peroxidases can also couple agmatine conjugates to form 
antifungal compounds known as hordatines in barley and 
(to a lesser extent) in wheat [59]. Finally, the blufensin-1 
gene was also induced in the mutant. This is a negative 
regulator conferring sensitivity to powdery mildew in bar- 
ley by modulating penetration resistance [60]. 

The repressed probe sets were related to systemic in- 
teractions with the environment, especially the ABA re- 
sponse, e.g. genes encoding dehydrins such as RAB 
proteins and other LEA proteins that accumulate during 
seed maturation and in vegetative tissues in response to 
abiotic stress [61]. RAB is a LEA family dehydrin that is 
induced by ABA, and the Arabidopsis ortholog RAB18 
is a useful marker for the loss of ABA sensitivity e.g. in 
plants that are deficient for the SWI/SNF chromatin- 
remodeling complex [62]. The gene is upregulated when 
a relevant repressor is mutated, e.g. SUMO E3 ligase 
SIZl [63]. We hypothesize that the mutant induces a re- 
pressor that acts against the ABA pathway or downregu- 
lates an inducer of the same pathway. Indeed, two 
transcription factors were repressed in the mutant under 
control conditions, the Q-type C2H2 zinc finger tran- 
scription factor TaZP23 and the homeobox-leucine zip- 
per protein ATHB-6. TaZP23 is upregulated in leaves 
and roots during drought stress and also following ABA 
treatment, suggesting a role in ABA-dependent gene regu- 
lation [64]. ATHB-6 may act as a growth regulator in re- 
sponse to water deficit, and is induced by ABA, salinity 
and dehydration stress [65]. It also regulates the ABA sig- 
naling pathway, acting downstream of ABIl [66]. 

The repressed probe sets also included two endonucle- 
ase genes related to the barley genes Bnuc and BENl, 
which are expressed during apoptosis [67,68] . Bnuc is also 
induced by salinity stress [69]. The papain-like cysteine 
proteinase (HvPAP14) and TdS40 were also represented 
in the repressed probe sets, and these also facilitate senes- 
cence and apoptosis [70]. HvPAP14 also generates mature 
forms of storage proteins in seeds, provides free amino 
acids during germination [71], and confers stress tolerance 



and defense against pathogens [72,73]. TdS40 is induced 
during post- flowering senescence in wheat [74], and also 
by JA, SA and in chlorotic/necrotic leaf tissue following 
infection with Pyrenophora teres [75]. 

There were few differences between the mutant and 
wild-type plants in the absence of water stress at anthe- 
sis and no statistically significant differentially-expressed 
genes could be identified also at the tillering stage (data 
not shown). Therefore despite a clear difference in 
chlorophyll content and chloroplast structure, the muta- 
tion appears to have a minimal impact on metabolism 
and gene regulation, a finding that does not provide suf- 
ficient explanation for the yield loss detected under field 
conditions. 

Water stress enhances the differences between mutant 
and wild-type plants 

We investigated the molecular basis of the difference in 
yield between the two genotypes by exposing them to 
water stress, which is the most common environmental 
constraint for durum wheat [4]. Most of the genes modu- 
lated by water stress were similarly regulated in both ge- 
notypes, suggesting a common response (Additional file 3: 
Table S2). Genes involved in the synthesis of osmopro- 
tectants were upregulated in both genotypes, whereas 
genes involved in ABA-dependent signaling, ethylene 
and JA synthesis were downregulated in both geno- 
types. Genes related to primary metabolism (amino acid 
and carbohydrate catabolism) were also similarly mod- 
ulated, indicating that both genotypes use metabolic 
strategies to counter the effects of drought. In contrast, 
genes encoding photosystem components and enzymes 
representing carbohydrate metabolism and the tricarb- 
oxylic acid (TCA) cycle showed significant differences 
that may explain the difference in performance between 
wild-type and mutant plants in the field. 

Photosynthesis and pigment synthesis 

Both transcriptomic and proteomic analysis indicated 
that the components of photosystems I and II (PSI and 
PSII) and ATPase were differentially expressed in the 
wild-type and mutant plants under water stress. Several 
PSI proteins were underrepresented in the mutant, in- 
cluding reaction center (RC) protein subunit IV (PSI-E), 
which is required for optimal electron transport to ferre- 
doxin and flavodoxin [76], and light harvesting complex 
LHCAl (CAB lB-21) (Table 9). PSI takes longer than 
PSII to recover from photo-inhibition because RC pro- 
teins are degraded instead of being repaired [77] . There- 
fore, the loss of PSI RC proteins in the mutant is likely 
to exacerbate the damage caused by stress. 

PSBO and PSBP maintain the structure of PSII [78] 
and were specifically overrepresented in the mutant. 
CP29 was upregulated in both genotypes, but this was 
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underrepresented in the mutant under stress conditions 
(Table 9). This protein binds LHCB proteins to the reac- 
tion center and confers photoprotection [79]. The LHCB 
proteins were downregulated in both genotypes and un- 
derrepresented in the mutant. In contrast, the genes en- 
coding LHCB proteins were upregulated in the mutant 
(Table 7). The beta and delta subunits of ATPase syn- 
thase were also overrepresented in the mutant, suggest- 
ing that the H* flux across the thylakoid membrane is 
greater in the mutant, with a knock-on effect on ATP 
production. 

Three genes responsible for the synthesis of carotenoid 
precursors were upregulated by stress specifically in 
the mutant: l-deoxy-D-xylulose-5-phosphate synthase 
(DXS), hydroxymethylglutaryl-CoA synthase (MVAl) 
and phytoene synthase (PSY). One gene encoding 
glutamyl-tRNA reductase 1 (HEMAl), which is respon- 
sible for the synthesis of chlorophyll precursors, was also 
upregulated by stress specifically in the mutant. The pig- 
ment content of the mutant plants did not increase rela- 
tive to the control but the activation of these genes 
could help to prevent a steeper decline in pigment con- 
tent under water stress conditions (Table 3). 

Redox changes in the chloroplast result in the stoi- 
chiometric adjustment of photosystem proteins and 
regulate the expression of both nuclear and plastid 
photosynthesis-related genes [80-82]. Because the mu- 
tant showed evidence of oxidative stress even under 
control conditions, the re-organization of the photosyn- 
thetic apparatus under water stress conditions could re- 
flect the additional oxidative stress that occurs during 
the early stages of drought. The loss of photosynthetic 
proteins (especially LHCB) by the mutant in response 
to drought stress reduces the amount of captured en- 
ergy and thus avoids the over-excitation of PSII which 
would generate additional ROS and cause oxidative 
damage [83]. 

Oxidative stress and ROS signaling 

The balance between ROS synthesis and scavenging de- 
termines the level of oxidative stress in plants [54,84]. 
ROS scavenging and detoxifying proteins such as peroxi- 
dases, glutathione transferases and glyoxalase were gen- 
erally modulated in a similar manner in both genotypes, 
but more probe sets of this group were downregulated 
specifically in the mutant under stress. In the same dir- 
ection, a glutaredoxin gene was upregulated in the wild- 
type plants under stress but was unaffected in the 
mutant (Table 7). Thylakoid ascorbate peroxidase and 
catalase were upregulated specifically in the mutant, sug- 
gesting a greater need for protection against ROS de- 
rived from chloroplast and peroxisome activity (Tables 7 
and 9). H2O2 can spread from the chloroplast to other 
organelles by diffusion through aquaporin-like PIP2 



proteins [85] which were downregulated in the mutant 
to a lesser extent than in wild-type plants (Additional file 
3: Table S2). These results indicate that the mutant ex- 
periences greater oxidative stress than wild-type plants 
under water stress conditions. 

Two probe sets related to ROS and stomatal closure 
were specifically modulated in the mutant under water 
stress. The specific downregulation of glyoxalase gene 
suggests that the mutant is less able than wild-type 
plants to detoxify methylglyoxal using reduced glutathi- 
one, thus allowing the production of ROS and stomatal 
closure [86]. Similarly, the specific upregulation of re- 
spiratory burst oxidase homologue D in the mutant 
(Table 7) suggests that the ABA-induced production of 
ROS in guard cells is promoted to encourage stomatal 
closure [87]. However, the differential expression of 
genes that regulate the stomatal aperture did not result 
in any differences between the genotypes in terms of sto- 
matal conductance under water stress. 

ROS can induce Ca^^-dependent signaling that triggers 
interactions between different calcineurin B-like (CBL) 
sensors and their CBL-interacting protein kinase (CIPK) 
targets [88]. Genotype-dependent differences were found 
in expression of these genes: whereas CIPK3 and CIPKIO 
were upregulated by stress in both genotypes, CIPK2 and 
CIPK6 were upregulated specifically in wild-type plants 
and CBLl and CIPK23 were downregulated specifically in 
the mutant. These two proteins cooperate with CBL9 to 
activate K*-channels and regulate stomatal responses [89]. 
CIPKIO and CIPK6 belong to the SnRK3-type family and 
mediate nutrient sensing, stress responses and ABA sig- 
naling [90,91]. As a result of altered expression of CBLs 
and CIPKs different Ca^*- dependent downstream re- 
sponses under water stress are expected between mutant 
and wild-type plants. 

MPK3, which was specifically downregulated in the 
mutant, is also involved in ROS signaling [25]. Once ac- 
tivated, it phosphorylates other stress-related proteins 
such as the transcription factors ZATIO and AZF2, and 
enzymes in the ethylene synthesis pathway such as ACS 
[92], all of which were downregulated by stress in both 
genotypes. ZATIO (STZ) is a transcriptional repressor 
responsive to chitin, which mediates the response to 
photo-oxidative stress and ROS [93]. 

Despite not showing damaged leaves or altered stoma- 
tal response, the data suggest that the pale-green mutant 
differs significantly from wild-type plants in terms of 
ROS sensing and signaling under water stress. ROS gen- 
erated in the chloroplast could move more easily inside 
the cells and elicit a signal triggering specific responses. 

Carbotiydrate metabolism 

Transcriptomic profiling revealed differences between the 
genotypes in the regulation of carbohydrate metabolism 
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under water stress conditions (Table 7). The synthesis of 
osmoprotectants such as rafFinose family oligosaccharides 
(RFOs) and trehalose was specifically upregulated in the 
wild-type plants, including the genes encoding galactoki- 
nase and raffinose synthases. RFOs protect membranes and 
scavenge ROS [94]. Trehalose-phosphate synthase (TPS) 
was also specifically induced in wild-type plants and thus 
repressed in the mutant under stress. Trehalose and 
trehalose-6-phosphate act not only as osmolytes to stabilize 
proteins and membranes, but also as signaling molecules to 
regulate metabolism. 

Two of the probe sets upregulated in the mutant 
plants encoded sucrose: sucrose 1-fructosyltransferase 
and sucrose:fructan 6-fructosyltransferase (SSTl and 
SFT6), which are responsible for synthesizing fructans 
from sucrose in the vacuole [95,96]. The mutant plants 
therefore appeared to promote fructan synthesis for re- 
mobilization as free carbohydrates, which are directly 
correlated with the water-soluble carbohydrate content 
in wheat [96]. Conversely, probe sets encoding cellulose 
synthase and cellulose synthase-like enzymes were up- 
regulated specifically in the wild-type plants and xylosyl- 
transferases were repressed in the mutant plants under 
stress (Table 7). These enzymes divert carbon to the syn- 
thesis of cell wall polysaccharides (cellulose, hemicellu- 
lose and xyloglucans) and their loss reduces growth and 
also inhibits photosystem activity and carbohydrate 
metabolism genes [97,98]. This shows the importance of 
cell wall structure not only for growth and turgor 
changes under stress but also for carbohydrate metabol- 
ism. Furthermore, sucrose, RFOs and fructans can dir- 
ectly or indirectly scavenge ROS, as reviewed by Keunen 
etal. [99]. 

The general picture that emerges from the modulation 
of carbohydrate metabolism suggests that the pale-green 
mutant is less able to deal with water stress than wild- 
type plants. Cell wall synthesis genes are activated in the 
wild-type plants whereas fructans are produced in the 
mutant. This correlates with the lower biomass of the 
mutant plants at anthesis and the lower yield under 
rain-fed conditions (Table 2). 

TCA cycle and related pathways 

The tricarboxylic acid (TCA) cycle generates energy 
and reducing power (ATP and NADPH) but consumes 
H2O and releases CO2. The substrates feeding the cycle 
are acetyl-CoA and other intermediates from lipid or 
protein catabolism. Acetyl-CoA is also produced from 
pyruvate (derived from glycolysis) by the action of the 
pyruvate dehydrogenase complex, and the subunit El 
gene representing this complex was upregulated in both 
genotypes under water stress. However, under condi- 
tions of carbon starvation or demand, additional strat- 
egies are used to supply TCA substrates [100,101]. The 



wild-type plants appear to be more effective at these 
strategies given the specific upregulation of the follow- 
ing genes (Table 7): 

i) Genes encoding products that convert amino acids 
into substrates suitable for the TCA cycle [102-104] 
were upregulated in both genotypes (Additional file 
3: Table S2), but three were upregulated specifically 
in wild-type plants: methylcrotonoyl-CoA carboxyl- 
ase, a and p subunits (MCCA and MCCB), and 2- 
oxoisovalerate dehydrogenase E2 (BCE2). 

ii) Genes encoding the NADP-malic enzyme (NADP- 
ME 3 and NADP-ME 4), which converts malate into 
pyruvate, and phosphoenolpyruvate carboxykinase 
(PCKl), which converts oxaloacetate into phospho- 
enolpyruvate, were upregulated specifically in the 
wild-type plants. 

iii) Genes encoding isocitrate lyase (ICL) and malate 
synthase (MLS) from the glyoxylate cycle were 
upregulated specifically in wild-type plants and thus, 
repressed in the mutant. This cycle produces succin- 
ate and malate acting as a bypass to avoid the loss of 
CO2 from the normal TCA cycle. It also metabolizes 
glyoxylate generated by photorespiration and pre- 
vents photo-inhibition [105]. 

iv) Genes encoding glutamate dehydrogenase and 
glutamate decarboxylase (GDH and GDC) were 
specifically downregulated and consequently 
repressed in the mutant under stress. They 
participate in glutamate metabolism, which balances 
nitrogen and carbon use in the cell [106] either by 
storing nitrogen as GABA, glutamate or glutamine, 
or by supplying carbon back to the TCA cycle. A 
functional GABA shunt is required for stress 
tolerance [94]. In the pale-green mutant, carbon 
stored as glutamate and GABA may not be intro- 
duced efficiently into the TCA cycle, disrupting the 
nitrogen/carbon balance and reducing growth under 
carbon demand or starvation imposed by stress 
[107]. GDH can also be inhibited by excess light or 
sucrose [108]. 

The feed-in to the TCA cycle and the ability to derive 
energy from it appear to be more limited in the mutant 
under stress, potentially explaining the differences be- 
tween genotypes in terms of biomass and yield in the 
field experiments. The mutant mobilizes carbon less effi- 
ciently, with a significant impact on grain filling and 
yield as observed in our field experiments. 

Regulators of carbon metabolism and signaling 

Genes that regulate carbohydrate and nitrogen metabol- 
ism were also differentially expressed between the wild- 
type and mutant plants under water stress, helping to 
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explain the basis of the metabolic differences between 
the genotypes (Table 7). 

The BZIP63 gene was specifically downregulated and re- 
pressed in the mutant plants under stress. The transcription 
factor encoded by this gene is an important node in the 
glucose-ABA interaction network and is repressed by both 
ABA and glucose [109]. It may promote starvation toler- 
ance by regulating transcription of genes such as the raffi- 
nose synthase (RS6, upregulated in the wild-type plants). 

AP2/EREBP proteins mediate hormone, sugar and 
redox signaling in the context of cold and drought stress 
[110]. Each member has a unique expression profile and 
the individual proteins are differentially regulated by glu- 
cose [111]. We found that RAP2.12 and RAP2.7 were 
specifically upregulated in wild-type plants and as a re- 
sult appeared repressed in the mutant, and that RAP2.8 
was specifically downregulated and also appeared re- 
pressed in the mutant, suggesting that the mutant and 
wild-type plants respond differently to the presence of 
sugars. 

PI3P5K was upregulated specifically in the mutant plants. 
This enzyme increases the levels of phosphatidylinositol-5- 
phosphate, which functions as a signaling molecule during 
osmotic stress and can repress the expression of WRKY70 
[112]. SALl, which was downregulated specifically in the 
mutant, acts as negative regulator of drought tolerance by 
modulating the levels of different sugars and controlling 
the corresponding signaling pathways [113]. 

Since the results highlighted some modifications of sig- 
naling molecules, the enzymes of the metabolism could be 
fme-tuned specifically in each genotype to adjust the con- 
centration of metabolites. 

Conclusions 

Under well-watered glasshouse conditions the reduced 
chlorophyll and carotenoids content of our pale-green 
mutant resulted on mild oxidative stress symptoms that 
could be compensated without altering photosynthesis 
and physiological performance by the induction of stress- 
response genes. When the availability of water was limited, 
we observed more extreme modulation of antioxidant en- 
zymes and photosynthetic proteins in the mutant than in 
wild-type plants, suggesting that the latter were better 
protected against oxidative stress. Concomitant molecular 
changes in the mutant revealed that the regulation of 
carbohydrate metabolism differed from wild-type plants, 
particularly the energy-producing pathways and pathways 
feeding the TCA cycle. The modification of carbohydrate 
metabolism in the mutant resulted in lower biomass accu- 
mulation and hence reduced yield under field conditions. 
This yield reduction was significant (10%) under irrigation, 
but was exacerbated in rain-fed conditions (14%) in agree- 
ment with the gene expression study, thus indicating that 
the mutant plants were more sensitive to water scarcity. 



Our results support the conclusion that the pale-green 
mutant was less able to adapt to terminal water stress, 
sustaining that reduced pigment content may be disadvan- 
tageous in durum wheat under water limited conditions, 
although a case by case study would be required when 
working with other pale-green mutants. 

Methods 

Plant material 

We compared two durum wheat {Triticum turgidum L. 
var. durum) genotypes: the wild-type variety Borgia 
(from the cross IRTA-1004 x Bidi 17) and a derived 
pale-green mutant (MD-597). The mutant was selected 
in 2002 in the M2 generation after treating the wild-type 
variety with the mutagen sodium azide (NaNa) and back- 
crossing twice with the parent to remove other muta- 
tions. Fi offspring of backcrosses with the wild-type 
were similar to it, while the ratio of the character segre- 
gation was 3:1. Therefore, the mutation is most likely 
controlled by a recessive nuclear gene.The mutant plants 
used in our experiments were the stable eighth gener- 
ation individuals (Mg) derived from the isolated original 
mutant plant (Figure 4). 

Experimental treatments 
Glasshouse experiment 

We used a split plot design with three replicates (eight 
pots per replicate) under controlled conditions in a 
glasshouse. The main plots were the water treatments 
(well-watered and water-stressed) and the subplots were 
the genotypes (wild-type and mutant). Two vernalized 
seedlings of the same genotype were transplanted into a 
4.0-L pots (19 X 15 cm) filled with a blend of white and 
frozen through black sphagnum peat mixture (Klas- 
mann-Deilmann GmbH), pH 6.5, conductivity 40 ms/m, 
NPK 14:16:18. Each genotype was planted in 48 pots, 
giving a total of 96 plants. All plants were grown with 
the following maximum/minimum temperatures and 
photoperiods: 13°C/7°C, 12 h daylight from emergence 
to the first node detectable stage, 18°C/10°C, 14 h day- 
light from then to heading, 20°C/15°C, 16 h daylight 
from heading to anthesis, and 25°C/20°C, 18 h daylight 
after anthesis. The pots were rotated twice a week to 
avoid temperature and light position effects. 

Both genotypes were grown at two levels of soil water 
content: i) control, in which the soil was maintained 
at field capacity during the whole experiment, and ii) 
terminal water stress, in which irrigation was field cap- 
acity from sowing to booting (Zadoks stage 45; [114]) 
and then reduced until maturity by watering the pots 
every second day with the amount of water necessary to 
maintain the pots at 1/3 of field capacity. The stress 
treatment was initiated by withholding water at the an- 
thesis stage. The soil water content was monitored every 
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second day using soil moisture probes (ECH2O dielectric 
aquameter, Decagon Devices, Inc. Pullman, WA, USA) 
in the center of the pots and was maintained by weigh- 
ing pots daily and watering them accordingly. Data of 
the probes indicated that field capacity conditions corre- 
sponded to a soil water content around 50 ml water/ 
100 ml of freshweight soil, while the water stress treat- 
ment reduced soil water content until ~13 ml water/ 
100 ml soil (Table 1). The flag leaves of four plants per 
experimental unit were measured between 11 a.m. and 
1 p.m. solar time at the booting and anthesis stages to 
determine (Table 1): i) the midday leaf water potential 
(Tw) using a pressure chamber (Model 3005; Soil Mois- 
ture Equipment, Santa Barbara, CA); ii) leaf conductance 
(gs) using a steady-state porometer (LI-1600, Li-Cor, 
Lincoln, NE); and iii) net CO2 assimilation rate (An) 
using a portable infrared gas analyzer (Model LCA-4: 
Analytical Development, Hoddesdon, UK). For the ex- 
pression profiling experiments, the flag leaf from five 
plants was collected for each of three biological repli- 
cates and for each water regime, between 11 a.m. and 
1 p.m. solar time. The tissues were immediately frozen 
in liquid nitrogen and stored at -80°C. Each flag leaf 
was ground in liquid nitrogen and 100-mg and 40-mg 
aliquots were prepared for transcriptomic and proteomic 
analysis, respectively. Before RNA and protein extrac- 
tion, five samples were pooled for each biological 
replicate. 

Field experiments 

Wild-type and mutant plants were grown in the field for 
three crop seasons (2010-2012) under irrigated (control) 
and rain-fed conditions in north-eastern Spain (latitude 41° 
38'N, longitude 0°23' E, elevation 250 m a.s.l.) Soil type was 
Mesic Calcixerolic Xerochrept with a silty-clay texture and 




Figure 4 Field plots of the pale-green mutant (IVlD-597) and the 
wild-type cv. Borgia. 



pH of 8.1. Plots were sown to a density of 500 seeds/m^ in 
randomized complete block designs with three replicates, 
each comprising eight rows (12.5 m in length and 15 cm 
apart). Field conditions are in Table 10. Irrigation was pro- 
vided four times by flooding the first crop season (on Dec. 
8th and at montiy intervals from March to May) and by 
sprinkler irrigations the successive years in which frequen- 
cies were determined according to the environmental con- 
ditions (Additional file 8: Figure S3). Fertilizer was applied 
as appropriate, and weeds, diseases and pests were con- 
trolled according to standard agricultural practices. Plant 
biomass was determined at anthesis using a sample of 
0.5 m row length and crop dry weight (g/m ) following the 
methodology described by Villegas et al. [115]. Plots were 
mechanically harvested at ripening and the yield was pre- 
sented at 12% moisture level. 

Microscopy 

Three plantlets at the three fully-developed leaves stage 
(Zadoks stage 13; [114]) from the glasshouse experi- 
ment, irrigated at field capacity, were collected from 
each genotype. Small segments (1 mm) from the central 
part of the third-leaf blades were fixed in 2.5% glutaral- 
dehyde in 0.1 M phosphate buffer (pH 7.2) at 4°C. After 
three washes in 0.1 M phosphate buffer (pH 7.4) the 
samples were incubated in 1% osmium tetroxide and 
2.5% potassium ferrocyanide for 2 h at 4°C. After ten 
washes for 10 min each with 0.1 M sodium acetate, and 
incubation for 30 min in 0.5% uranyl acetate followed by 
two washes for 10 min in 0.1 M sodium acetate at 4°C, 
the samples were dehydrated in increasing concentra- 
tions of ethanol followed by two washes with sodium 
acetate and finally embedded in epoxy EMBED-812 resin 
(Electron Microscopy Sciences) for three days at 60°C. 
Ultrathin sections (70 nm) contrasted with Reynolds' 
lead citrate stain were examined using a Zeiss TEM 910 
transmission electron microscope. 

Pigment measurements 

Pigments were quantified in four wild-type and mutant 
plants from each treatment regime harvested at anthesis. 
The blade of the flag leaf was dissected, weighed and imme- 
diately frozen in liquid nitrogen. The pigments from 
100 mg lyophUized tissue were extracted with acetone for 
20 min at 40°C, filtered and mixed with 9:1 hexane:ethyl 
ether and 2% NaCl. The organic fraction was evaporated in 
a nitrogen stream and argon was used for storage at -80°C. 
The pigments were separated on an ACQUITY Ultra Per- 
formance LC~ system (Waters, Milford, MA, USA) linked 
to a PDA 2996 detector (Waters, Milford, MA, USA) fol- 
lowing the protocol of Taylor et al. [116] but using a modi- 
fied gradient. 

Carotenoids were isolated using a YMC C30 carotenoid 
3-nm, 2.0 X 100 mm HPLC column (Waters, Milford, 
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Table 10 Details on field experiments 



Crop 


Water 


Sowing 


H3rv6StinQ 


Water input 


M63n t6mpGr3tur6 (°C) 


M6dn 3ir r6l3tiv6 humidity {%) 


season 


regime 


date 


date 


(rainfall -i- irrigation, mm) 






2009-2010 


Rainfed 


1 3 Nov. 2009 


12 July 2010 


287 


114 


70.0 




rrigated 


1 3 Nov. 2009 


13 July 2010 


487" 


114 


69.9 


2010-2011 


Rainfed 


26 Nov. 2010 


5 July 2011 


245 


11.5 


70.3 




rrigated 


27 Nov. 2010 


7 July 2011 


395^ 


11.7 


70.2 


2011-2012 


Rainfed 


9 Dec. 201 1 


2 July 2012 


160 


11.5 


62.8 




rrigated 


8 Dec. 201 1 


12 July 2012 


335b 


12.0 


62.9 



^Flooding irrigation. 
^Sprinkler irrigation. 



MA) with a mobile phase comprising solvent A (metha- 
nolrwater, 97:3, v/v) and solvent B (100% methyl tert-huty- 
lether), both solvents containing 0.05% triethylamine. The 
sample was warmed to 25°C and a volume of 10 |iL was 
injected into the column, which also was thermostatically 
maintained at 25°C. The flow rate was set to 0.5 mL/min. 
The gradient program was set as follows: initial conditions 
98% solvent A and 2% solvent B for 13 min; changed with 
a linear gradient to 62% solvent A and 38% solvent B in 
1 min; hold for 2 min; changed with a linear gradient to 
32% solvent A and 68% solvent B in 1 min; hold for 5 min; 
returned to initial conditions in 4 min, followed by equili- 
bration for 5 min. The data analysis was done with Mas- 
sLynx™ software version 4.1 (Waters, Milford, MA, USA). 
Carotenoids were identified by monitoring the order of 
elution from the column, the ultraviolet and visible spec- 
tra and the spectral fine structure and comparing them 
with standards and previous references [117]. Standards 
trans-p-apo-8'carotenal, chlorophyll a and chlorophyll 
b were purchased from Sigma-Aldrich (USA), whereas 
lutein (L), zeaxanthin (Z), P-carotene (^-car), neoxanthin 
(N), violaxanthin (V) and antheraxanthin (A) were ac- 
quired from CaroteNature (Switzerland). The data were 
analyzed by ANOVA with the SAS-STAT package and 
the means were compared with a Tukey test at 5% 
probability level. 

RNA isolation and array hybridization 

Total RNA was extracted using Trizol reagent (Invitrogen, 
Carlsbad, CA, USA) and cleaned with RNeasy MinElute 
(QIAGEN) columns following the manufacturer's instruc- 
tions. Quality was determined using an Agilent Bioanalyzer 
2100. The RNA samples were processed following the Affy- 
metrbc Genechip Expression Analysis Technical Manual. 
Single-stranded and double-stranded cDNAs were synthe- 
sized using the Affymetrbc GeneChip" 3' IVT Express 
Kit. Biotin-conjugated nucleotides were incorporated into 
aRNA by in vitro transcription to generate cRNA which 
was purified to remove unincorporated NTPs, salts, en- 
zymes and inorganic phosphate before fragmentation and 
hybridization onto 3' expression arrays. Specifically, we 



used the Affymetrbc GeneChip' Wheat Genome Array. The 
arrays were washed and stained on an ASymetrbc Fluidics 
Station followed by scanning with a GeneChip Scanner 
3000. 

GeneChip® quality analysis 

The quality of microarray hybridizations was determined 
by checking the following parameters: i) the standard Affy- 
metrbc controls to evaluate labeling and hybridization (B2 
oligonucleotides, PolyA controls such as lys, phe, thr and 
dap, and hybridization controls such as BioB, BioC, BioD 
and Cre); ii) RNA degradation plot relative to probe signal 
intensities of the Actin and GAPDH control genes, using 
the R library "simpleaffy"; iii) percentage of "Present calls" 
on the basis of the MAS 5.0 algorithm; and iv) R^ linear 
correlation coefficients among biological. 

Data processing and analysis 

Hybridization data were processed as recommended by 
Aprile et al. [118]: i) normalization of the raw data by 
Robust Multi-Array Average (RMA) [119] using the R 
package Affymetrix library [120]; ii) MAS 5.0 algorithm 
on raw data in order to produce a detection call for each 
probe set "Present", "Marginal" or "Absent" (all probe 
sets that did not show a "Present" call in all repetitions 
of at least one sample were removed); iii) the filtered 
data were imported into the software Genespring GX7.3 
(Agilent) for the analysis of differentially-expressed 
genes using a Welch i-test with Benjamini and Hochberg 
false discovery rate correction for multiple tests [121]. 

The differences in gene expression were considered to 
be significant at P < 0.05 and the induction or repression 
ratio was >2-fold. Principal component analysis (PCA) 
was then used to validate the role of the genotype and to 
verify the effect of water stress between the wild-type 
and mutant plants (Figure 2). 

The terms induced or repressed were used when de- 
scribing the comparison on gene expression between 
both genotypes independently of the water treatment 
(mutation effect). And the up/downregulated terms were 
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used when comparing the two water treatments in the 
same genotype (stress effect). 

In order to identify overrepresented gene classes within 
selected groups of genes compared to the entire array, the 
Arabidopsis thaliana best BLASTX annotations of the 
wheat probe sets were used as input for the MIPS Func- 
tional Catalogue Database (FunCatDB) [122], which pro- 
vides a statistical survey of the functional distribution of a 
given set of genes. Only the functional categories with a 
cut-off of P < 0.005 were considered. BLAST searches were 
carried out using Plant Expression Database (PLEXdb) 
[123] where a complete list of additional annotation for all 
probe sets on the wheat genome array is available. The ori- 
ginal microarray data have been deposited and are access- 
ible in PLEXdb (http://www.plexdb.org) as experiment 
TA48: "Pale-green durum wheat mutant under terminal 
drought stress condition" and in NCBI's Gene Expression 
Omnibus (NCBI GEO; http://www.ncbi.nlm.nih.gov/geo) 
as experiment GSE45563. 

Validation of array data with qRT-PCR 

Quantitative real-time RT-PCR (qRT-PCR) was used to 
confirm the differential expression of the probe sets iden- 
tified by microarray analysis. The total RNA isolated for 
the array experiment was treated with DNA-free'" DNase 
Treatment & Removal Reagents (Ambion, Life technolo- 
gies) and cDNA was synthesized using the iScript cDNA 
synthesis kit (Bio-Rad, Hercules, CA, USA). PGR amplifi- 
cation was carried out in a 25 |iL final volume containing 
2x IQ SYBR Green Supermtx (BioRad, Hercules, GA, 
USA), 500 nM of each specific forward and reverse 
primers and 10 ng of cDNA as the template, using the 
Bio-Rad GFX96 sequence detector system. The amplifica- 
tion program comprised a denaturation step at 95 °G for 
3 min followed by 40 cycles of 95°G for 10 s, 60°G for 30 s 
and 72°G for 20 s. The amplification products were vali- 
dated by melt curve analysis using standard curves with 
an efficiency of 90-100%. Several repetitions were used to 
establish the efficiency of the standard curve for the in- 
ternal reference gene. Efficiency was auto-calculated using 
the Bio-Rad GFX Manager software [124]. Gt value and 
normalized expression were calculated using the AAG(t) 
method [125]. The no reverse transcriptase and no tem- 
plate negative controls were used to ensure absence of 
contamination and of nonspecific amplifications. Primer 
specificity was further confirmed by gel electrophoresis. 

Specific qRT-PGR were used to verify nine probe 
sets corresponding to the following genes: Chitinase 3 
(TaA£Ex;.128418.43.Sl_at), Blufensin 1 (TaAfik.26815.1.Sl_at), 
Germin-like protein (Ta.l69.1.Sl_x_at), Naringemn,2- 
oxoglutarate 3-dioxygenase (TaAf6c.6092.1.Sl_at), WRKY45- 
like transcription factor (Ta.8614.2.Sl_x_at), Agmatine cou- 
maroyltransferase (Ta.8228.1.Sl_at), RAB protein (Ta.2638. 
l.Sl_at), Q-type C2H2 zinc finger protein (ZFP23) (TaAfEx. 



98394.1.Sl_at) and Dehydrin {WZYl-1; Ta.l3255.1.Sl_at). 
Primers were designed near the 3', according to GenBank 
(NGBI) database sequences (Additional file 9: Table S6). 
Data were calculated from the calibration curve and 
normalized using the expression curve of the Actin tran- 
script (GQ339780) corresponding to probe set Ta.28253.1. 
Sl_at, selected among several probe sets with low vari- 
ation between samples as the internal reference gene 
(GV = 0.036). The qRT-PGR data were compared to the 
corresponding microarray expression values by mean 
of Pearson product-moment correlation coefficients 
(Additional file 1: Figure SI). 

Protein extraction 

Protein was extracted as described by Jiang et al. [126] with 
modifications. The powder was resuspended in 2 mL ice- 
cold 10% trichloroacetic acid in acetone containing 0.07% 
dithiothreitol (DTT) and protease inhibitor cocktail (P9599 
Sigma). After overnight incubation at -20°G, the mixture 
was centrifuged at 4 °G for 30 min at 14000 rpm. The pellet 
was washed four times with acetone containing 0.07% DTT 
and centrifuged for 5 min at 14000 rpm. The pellet was left 
to stand on ice until the acetone evaporated completely. 
Protein Extraction Reagent type 4 buffer (G0356 Sigma) 
was used to resuspend the pellet overnight at room 
temperature while mixing. Proteins were separated from in- 
soluble components by centrifugation for 20 min at 
14000 rpm and were stored at -80°G. The Bradford protein 
quantification method [127] was used to quantify the pro- 
teins using bovine serum albumin as a standard. 

2-D Electrophoresis and detection 

For the first dimension separation, samples were mixed 
with rehydration buffer (7 M urea, 2 M thiourea, 1% 
G7BzO detergent, 40 mM Trizmabase, 50 mM DTT, 1% 
IPG buffer pH 3-10, and 0.002% bromophenol blue) to a 
total volume of 315 |iL. Actively rehydrated (50 V/gel, 
20° G) 18-cm Bio-Rad IPG strips (pH 5-8) were loaded with 
protein and separated by isoelectric focusing using a Pro- 
tean lEF Gell system (Bio-Rad) following the manufacturer's 
instructions. The sequential gradient procedure was 500 V/ 
linear for 30 min, 1000 V/linear for 1 h, 10000 V/linear for 
1 h, and 10000 V/rapid until a total of 55,000 Vh. The 
current limit was 50 |iA per IPG strip. After lEF separation, 
the gel strips were incubated for 15 min in the equilibration 
buffer (375 mM Tris base, 6 M urea, 20% glycerol, 2% SDS) 
containing 2% DTT, followed by 15 min in the same buffer 
containing 2.5% iodoacetamide instead of DTT. For the 
second-dimension separation, SDS-PAGE with 11% 22 x 
20 cm polyacrylamide gels were loaded and separated in an 
Ettan DALTsix Electrophoresis Unit (GE Healthcare) with 
0.25 M Tris-HGl (pH 8.8), 1.92 M glycine, 1% w/v SDS 
electrophoresis buffer, at 8 mA/gel overnight. Three repli- 
cate analytical gels, with 50 |ig of proteins each, were run 
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for each protein extract. The preparative gels were loaded 
with 350 i^g of protein. After electrophoresis, the analytical 
and preparative SDS-PAGE gels were fluorescendy stained 
with Flamingo or Oriole (Bio-Rad) respectively, and 
scanned with a Versadoc MP 4000 system (Bio-Rad). 

Image analysis 

Gel imaging and spot quantification were carried out 
using the PD-Quest Advanced 2D Gel Analysis software 
vS.O.l (Bio-Rad) applying automatic spot detection and 
matching followed by manual/visual validation. The 
image was normalized by local regression and a scaling 
factor of 10 (ppm) was applied to avoid non-expression 
related variations in spot intensity. Spots showing more 
than 2-fold changes in intensity with a coefficient of 
variation <50% and with statistically significant differ- 
ences by Student's t-test (P < 0.05) were defined as differ- 
entially expressed. 

Protein digestion and identification by mass 
spectrometry (MS) 

The spots of interest were excised manually from pre- 
parative gels and digested with trypsin according to 
Perez-Hedo et al. [128]. We then applied 1 ^iL of tryptic 
peptide solution to a MALDI plate, dried it at room 
temperature and covered it with 1 |iL of saturated a- 
cyano-4-hydroxycinnamic acid prepared in 50% v/v 
ACN containing 0.1% TFA. 

Proteins were identified using an AutoflexSpeed 
MALDI-TOF/TOF mass spectrometer (Bruker Daltonics) 
by peptide mass fingerprinting (PMF). For those spots lack- 
ing identification, 3-4 peptides were analyzed by MS/MS 
mode. Mass spectra (mode reflectron, MH^) were acquired 
by FlexControl v3.0 software (Bruker Daltonics), recording 
in the range 800-4500 Da, and the MS/MS information 
was obtained in LIFT (Laser-Induced Forward Transfer) 
mode. MS spectra were externally calibrated using Peptide 
Calibration Standard II (Bruker Daltonics). The peak lists 
obtained were compared against the Swiss -Prot, Trembl 
and non-redundant NCBI protein databases, using the 
MASCOT software package v2.3 (Matrix Sciences, UK). 
The search parameters were set as monoisotopic peptide 
masses, carbamidomethylation of cysteine and oxidation 
of methionine as fixed and variable modifications, respect- 
ively, one trypsin missed cleavage and a maximum of 
±100 ppm for PMF peptide tolerance and ±0.4 Da for MS/ 
MS tolerance. Protein identification was accepted for those 
results showing a probability-based MOWSE score with 
P<0.05. We also analyzed MS and MS/MS combined 
spectra prepared using BioToolsv3.1 (Bruker Daltonics). 

Availability of supporting data 

The data set supporting the results of this article is available 
in the PLEXdb repository as experiment TA48. 
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